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ABSTRACT 


For many astrophysical situations, such as in solar flares or cosmic gamma- 
ray bursts, continuum' gamma rays with energies up to hundreds of MeV have been 
observed, and can be interpreted to be due to bremsstrahlung radiation by relativistic 
electrons. The region of acceleration for these particles is not necessarily the same 
as the region in which the radiation is produced, and the effects of the transport of 
the electrons must be included in the general problem. Hence it is necessary to solve 
the kinetic equation for relativistic electrons, including all the interactions and loss 
mechanisms relevant at such energies. The resulting kinetic equation for non-thermal 
electrons, including the effects of Coulomb collisions and losses due to synchrotron 
emission, has been solved analytically in some simple limiting cases, and numerically 
for the general cases including constant and varying background plasma density and 
magnetic field. New approximate analytic solutions are presented for collision dominated 
cases, for small pitch angles and all energies, synchrotron dominated ^ases^both steady- 
state and time dependent, for all pitch angles and energies, and for cases when both 
synchrotron and collisional energy losses are important, but for relativistic electrons. 
These analytic solutions axe compared to the full numerical results in the proper limits. 
These results will be useful for calculation of spectra and angular distribution of the 
radiation (X-rays, 7-rays and microwaves) emitted via synchrotron or bremsstrah- 
lung processes by the electrons. We shall examine these properties and their relevance 
to observations in subsequent papers. 

Subject headings: plasma-radiation mechanisms-Sun:flares-Sun:X-rays-X-rays:bursts 
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I. Introduction 


In many astrophysics! situations the observed electromagnetic radiation is produced 
by accelerated electrons with non-thermal or non-maxwellian distributions (typically 
with power law energy spectra and anisotropic momentum distributions). Interaction 
of these electrons with ambient plasma generally with varied particle, photon and 
magnetic field densities produces the observed radiation through synchrortron, compton 
or bremsstrahlung processes. In general the acceleration site, normally a region of 
low density, high plasma turbulence or electric field, can be different than the region 
where the bulk of the radiation is produced. The distribution of particles derived from 
modeling of the emission process is not necessarily that of the accelerated particles but 
is modified during the transport from one region to another. The two distributions are 
related by the particle kinetic equation. It is imperative then to understand the effects 
of transport in order to determine the distribution of the accelerated particles and to 
gain insight into the acceleration mechanism. 

In most analyses of non-thermal sources the transport effects are either ignored 
or treated in an approximate manner primarily because of the complexity of the 
problem. Many interactions with the ambient plasma such sis Coulomb collisions, 
inverse Compton scattering, Synchrotron cooling and interaction with both small scale 
electromagnetic field fluctuations (plasma waves) and large scale electric and magnetic 
fields can be simultaneously be important. Some simplfied cases have been analyzed. 

For example, in recent studies the synchrotron effects were considered by Lamb and 
Brainerd (1987), Coulomb collision effects by Leach and Petrosian (1981, hereafter LP) 
and Petrosian (1983), and synchrotron and magnetic field effects by Ho (1986) and 
MacKinnon and Brown (1988), In this paper we consider the effects of all three of these 
processes. In subsequent papers we shall apply the results from this study to microwave, 
X-ray and 7 -ray production in solar flares, and X-ray to 7 -ray production in gamma ray 
bursts. 

We use the Fokker-Planck method for solution of the kinetic equation and 
determination of the distribution f(x,p,t) in phase space. Because of the presence of 
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strong fields in general the particle drift across field lines is neglegible so that only the 
spatial coordinate s, the distance along the field line, and two components of momentum 
(parallel and perpendicular to the field lines) are needed. In most cases discussed 
below we replace the two momenta with the kinetic energy E and fi = cos a where 
a is the particle pitch angle. We shall also consider the steady state case which is a 
good approximation when the time scale of the modulation of the observed radiation 
is longer than the typical time scales of the transport processes. Thus we solve the 
kinetic equation for the distribution /(£J, where / dEdfids is the number density 
of particles. We utilize the Fokker-Planck treatment of the collisional effects and neglect 
the inverse Compton effects (which are in many ways similar to synchrotron effects), the 
effects of electric fields, plasma turbulence and self- absorption of synchrotron radiation. 

The fully relativistic equation including collisional and synchrotron effects as well as 
the effects of the field inhomogeneities can then be written as 
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where $ = ///?, 7 = E + 1 is the total energy, /9c = cyj\ - I/7 2 is the electron velocity, 
B is the magnetic field strength and E is a source term for the injected electrons. 

[See, e.g., LP, Leach (1984), Petrosian (1985) and McTieman (1989). Note that the 
convergence term is different here than in the equations in LP and Leach (1984), with 
the factor of (1 - y?) being inside the derivative d/dfi. This results in differences of less 
than 5% in comparison with the earlier results.] 

The steps leading to equation (1) and the definition of the collisional energy loss 
and diffusion coefficients C and C' = £C, and the synchrotron coefficient 5 are given in 
LP and Petrosian (1985). For a background plasma of fully ionized hydrogen, £ = 1 and 

C=4^nh,A = 2xlO-’(!|^)( I5; ^p) cm - 1 . (2 


where r 0 = e 2 /m e c 2 is the classical electron radius, n is the ambient proton or 
electron density, and In A is the coulomb logarithm. This simple relation is not true for 
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a partially ionized plasma or a neutral gas. In these cases the ratio £ depends on the 
energy; e.g., for a neutral background £ varies from ~ 1/12 to ~ 1/8 for energies from 


10 keV to 10 MeV. 

The synchrotron energy loss and pitch angle change terms axe proportional to 
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The coefficients S, C, and C' have the units of inverse length and are useful scales. We 
will find it convenient to define the ratios R 0 S S/C (f? 0 = 3230 B 2 /n for In A = 20), 
and 

R = /?V(1 - S)Ro • (4) 


as a measure of relative importance of the synchrotron and collisional energy losses (see 
Table 1). 

For the source term £ we shall assume that the electrons are injected at one point, 
the origin of the spatial coordinate s = 0 , so that E oc £(s) is zero everywhere 
except at s = 0. Consequently we solve the equation with E = 0 and use the injected 
distribution as a boundary condition. Furthermore, we will present the results in terms 
of the particle flux integrated over the cross-sectional area of the loop A, given by 
F = 0c f A = 0 2 c$A, which has units of s _1 keV -1 ster -1 . The source, or the boundary 
condition at s = 0, is assumed to have the form 


F(E, = 0) = F 0 (E)G(fi) = 2a- 2 F o (£0e-<“’-"?>/ a ’ . (5) 


Here a 2 is the dispersion in pitch angle and c*i is the peak direction. For lower values 
of a\ and a, the distribution becomes beamed along the field lines. For aj = 7t/2 and 
small a 0 the distribution wiU be of pancake form. In some cases we shall replace a 2 with 
sin 2 a = (1 — n 2 ), and when necessary we shall assume a power law energy spectrum 
given by 

F 0 (E) = KE~ S . (6) 

In the next section we discuss some analytic solutions for equation (1). In Section 
III we describe the results from numerical solutions of this equation for a variety of 
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injected pitch angle distributions and field configurations. A brief summary is presented 
in Section IV. 


II. Analytic solutions 

In some limiting cases, we find that analytic solutions are possible. The analytic 
solutions are useful in many ways. They can give good quantitative estimates within 
their domain of applicability; they can be used to test the accuracy of the complex 
n um erical results, and they provide a qualitative guide for more complicated cases 
outside the range of applicability. 


A. Collision dominated solutions (R <£. 1) 

For high densities and low values of the field strength and electron energy, the 
electron transport is dominated by Coulomb collisions and we can set S = 0. As shown 
by LP analytic solutions are then possible for small pitch angles. If a 2 1, the injected 
electrons are strongly beamed along the field lines and we can set n = (1 — a 2 / 2) and 
(1 _ n 2 ) = a 2 in equation (1). If we ignore the 0(a 2 ) terms, and define a normalized 
collisional column depth r c and energy parameter rj by 

dr c = Cds , dr) = fdE (7) 


then the flux at any depth is given by equation (18) of LP: 

The pitch angle distribution is gaussian at all r c with a dispersion a c given by 

a 2 =a 2 + C(£,r c ), 


( 8 ) 

( 9 ) 


where 


C(£,r c ) = 2£ln 


' S(q + r c ) 

. B(t 7 ) 
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From equation (7) we have r/ = E 2 /(E + 1) so that 
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Thus E(r]) = E and [/?(x)] 2 = 1 - l/(x + l) 2 relates the velocity 0 to energy. 

1) N on-relativistic limit. For non-relativistic particles E C 1 and r) as E 2 . For 
r c <C 1, which will be the case for these particles because they lose most of their energy 
by t c ~ rj ~ E*, 

E (>7 + r c ) = E^l+^JE 2 


and 

C(f:,r c ) = ein(l + r c /£; 2 ). (12) 


As shown by LP this solution is a good approximation up to very large values of r c / E 2 
and for injected pitch angle distributions with values of a 2 up to 0 . 40 ; much larger than 
expected considering the assumptions made. 


2) Relativistic limit. For the extreme relativistic (E > 1) electrons we have tj 
and E(tj + r c ) = E + t c . The diffusion in pitch angle is small and according to 
equations (9) and (10) the dispersion a 2 does not change appreciably with depth: 

4£t c 


a 2 r=« 2 o + 


E(E + t c ) ' 


(13) 


This implies an increase in dispersion with depth from a 2 to a 2 4- 4 £/E which is a small 
effect, except for highly beamed injection with a\E <C 1 . 

Equation ( 13 ), however, often overestimates the dispersion in pitch angle. We 
obtained the approximate solution equation (8) by setting /i = 1 in front of the d$j ds 
term on the left hand side of equation ( 1 ). This is reasonable for the non-relativistic 


case where the neglected term is of order a 2 and is insignificant in comparison with the 
effects of the diffusion term which is of order unity. As the electron energy increases the 
diffusion term becomes smaller and for sufficiently high energies (FJa 2 > 1) it becomes 
comparable to or smaller than the 0(a 2 ) term neglected. Thus for relativistic electrons 
we need a more accurate treatment of the n(d$/dr c ) term. As shown in part 1 of 
Appendix A in this limit the diffusion term can be treated as a perturbation leading to 
the approximate solution (A.4) for extreme relativistic electrons and all angles, 

= {wgrhd) w-wjm ■ < 14) 
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and to a more accurate solution for intermediate and high energies but for the small 
pitch angle regime (eq. [A. 12]); 


F(E,h,t c ) 


( m 

Kmv + r, 


0[E] V 2e- a ’/*? rrv ,, 
TTTTu — - — W? + ’’.)] . 


where a 2 c is given by equations (9) and (10) and the dispersion as a function of depth is 


d? = 


a r 2 

C ^ 


c [l+ 6 e T c a 2 c / 2(2 + r 1 + T c )] <ac ' 


Here 8 e = d\nF 0 {E)/d\nE\( E +r c ) and is equal to the spectral index 8 if F 0 is a 
power law. Note that equations (15) and (16) reduce to the non-relativistic limits of 
equations (8) through (12) in the proper limit i) ~ r c < 1. Hence we may use this 
corrected solution for all energies. 

Table 1 gives the values for a 2 and al along with a 2 , the dispersion obtained from 
numerical solutions of equation (1) including only the collision terms, but without the 
small pitch angle approximation, for four values of a 2 at r c /»7 = 1. For a, = 0.04 the 
diffusion effects are more important (4 /Ea 2 « 5 and 8a 2 / 2 « 0.10). Consequently, 
the dispersion increases with depth, while for a* = 0.40 the reverse is true (4/ Ea] « 

0.5 and 8a 2 / 2 ss 1.0) and the dispersion decreases. In all cases except for isotropic 
injection (a% = 40), a 2 provides an excellent approximation to the dispersion a 2 . For 
al >> 1 it is obvious that a 2 is not a good approximation but a 2 c provides a reasonable 
approximation. Figure 1 shows a comparison of the pitch angle distribution from the 
above analytic approximation (eqs. [15] and [16]) with that from the exact numerical 
solution for a power law injected flux given by equations (5) and (6) with 8 = 5 and 
a 2 = 0.40. The agreement is excellent at small pitch angles, and reasonable at somewhat 
higher pitch angles. This comparison gives an indication of both the accuracy of our 
numerical code and of the usefulness of the small pitch angle approximation. 

8 ) Flux integrated over pitch angle. In certain problems, for example those with 
straight magnetic field lines, and processes with isotropic cross sections, knowledge of 
the pitch angle distribution is not necessary. We define total flux of electrons of a given 
energy at a given depth to be F^E, t c ) = J +l F(E, r c )dfi and integrate the solutions 
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from equation (15) to obtain 
Fp(E,r c ) = 


/PF.[E(n + T C )] 


(17) 


P 2 [E(rj + T c )][l + 6 e r c a 2 / 2(2 + rj + r c )] 

which, as shown by LP, is independent of a 2 for non-relativistic electrons for all values 
of aj. For relativistic electrons, F M does depend on a 2 due to the presence of a 2 in the 
correction to the dispersion in equation (16). 

The qualitative behavior of is similar for all energies. The flux is constant from 
r c = 0 to r c ~ r) and afterwards decreases with increasing depth. For a power law 
injected flux, F 0 (E) = KE ~ fi , we find 

f KE(E 2 + r c )-(^>/ 2 for E « 1 (18) 

^^\K(E + T c )- 6 /[l+ST c a 2 c /2(E + T c )} for E > 1 ’ V ' 

For large values of t c > 17, the flux falls off as r c _(tf+1)/2 for E < 1 and as r~ s for 
E > 1. 

4) Spatially integrated flux. For a magnetic field which is both uniform and straight 
a useful quantity is F r (E, //), the flux integrated over column depth, which is defined as 
F r (E, fi) = / o °° F(E, / 1 , T c )dr c . This quantity can be sufficient for the study of spatially 
unresolved sources. The zeroth order approximation, equation (14), integrated over r c 
gives 

F r (E, h ) = mGM n^WE' . (19) 

J E 

We cannot obtain an analytic expression for F r by integrating equations (8) or (15) 
due to the complex r c dependence of the dispersions a 2 and a 2 . It is much easier to 
integrate the original equation over 0 < r c < 00. The integral of the source term E is 
equal to 

f 00 Sdr c = cF 0 (E,fi) = cp 2 * 0 (E,v) (20) 

J n 


and the equation becomes 


/? 2 7 2 — 


d_ 

dE 




2, Mr 


dn J 


-Hfl* 7 2 *o = 0, 


( 21 ) 


where $ r {E, fl) = / o °° $(F, fi, T c )dr c . 
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For relativistic energies the pitch angle diffusion term is small and as above we can 
treat it as a perturbation. We expand $ r in terms oi\/E and include the first order 
correction due to diffusion. The zeroth order solution for the flux is (since F r = c$ r for 
£> 1) 

F r {E, ft) = n f X F 0 (E', »)dE' , (22) 

J E 

and the first order solution is 

F,(E,v)-tL J°° F.(E‘, h )dE' 


(23) 


For FJE, n) = FJE)G(ji) this gives 

F r (E,ri = ?G(f 0 J~ F.(E')dE' ( 


S T E^GMd^ 


.2- 9 




)' 


(24) 


where we have defined 


E /“ E'-'dE' FJE")dE" 


/, K(E’)dE' 

which is a measure of the energy spectral index (for a power law F 0 , S T =6). 

For the small pitch angle regime with the injected gaussian distribution given in 
equation (5), we find 


(25) 


F r (E,») = 


2e‘ 


+ 


' 2/ °* 1 l 

H<*’/2)J J E 


where 


a 2 = 

E 


at 


/ F 0 (E')dE' , 

(26) 

> r E ‘ 

(27) 


For the power law injected energy distribution S T = 6 and because the dependence 
on energy of a\ is small (E > 1), the spatially integrated flux should have a power law 
index equal to 6 — 1. 

5) The total energy spectrum. Integration of F M over r c or F r over \i gives us the 
total spectrum F TOT (E) relevant for situations with isotropic processes and spatially 
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unresolved observations. This solution is obtained by integrating equation (19) over dfi. 
Thus 

F tot (E) = f' MMn J\y?F 0 (E')dE' <x ( « \ , (28) 

where the last relation is applicable fo a power law injected spectrum. In Section III we 
will compare these analytic results with numerical ones. 


B. Synchrotron dominated case (R \) 

In the opposite limit of high magnetic fields and particle energies synchrotron losses 
dominate. A general time dependent solution for the case with a uniform magnetic 
field and including synchrotron losses ( dlnB/ds = C = 0) is given in part 2 of 
Appendix A. For the steady-state case with continuous injection the result is given by 
equation (A. 24). This solution is valid for electrons of all energies, but it takes a simple 
form for relativistic energies. The behavior of non-relativist ic electrons is qualitatively 
similar to that of relativistic electrons. 


1 ) Relativistic limit. The relativistic limit of equation (A. 24) is given by 
equation (A.30), which for an injected flux of the form F(E, /u,0) = F 0 (E)G{fj.) reduces 
to 


F(E,fi,s) = 


F 0 [E/(l-r M /r ser )]G(it) 


(1 - T a /T, cr ) 2 

where we have defined a dimensionless depth r, = sS and 


r(E,ll)~ ((j_^ 2 )) E 


(29) 


(30) 


At a given pitch angle, F -» 0 at r lcr which decreases with increasing energy so that 
higher energy particles are stripped from the beam at smaller depths. The depth s cr 
corresponding to r $cr also decreases with increasing magnetic field strength through the 
B 2 dependence of 5. 

At a given t s , the flux becomes zero at a critical pitch angle a cr = cos 1 {^cr) given 
by 

Hcr(E, T 9 ) = + 4r » 2£?2 “ l ) ' ( 31 ) 
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Note that n cr increases as r, increases, and approaches 1 asr, — ► oo. Electrons with 
higher pitch angles are stripped from the beam as the depth increases. 

In order to see the initial trend of the distribution it is instructive to consider the 
small pitch angle regime for r $ •< r scr = (a 2 .E) -1 . In this case, we let n = 1 — a 2 /2 and 
G(n) = 2a~ 2 e~“ J / or ® to obtain 

F(E, fi, s) = 2 a; 2 F 0 (E)e- a3 '°' , (32) 

where (unless S e < 2) 

a]=al[l + (6,-2)alT.E]- 1 , (33) 

with S e = din F 0 (E)/d In E. The dispersion a] decreases with increasing energy, depth 
and magnetic field. The effects of the pitch angle term and the first order correction 
to the extreme relativistic approximation will add terms of order 1/E in the square 
brackets. 


3 ) Spatially integrated flux. The general expression for the flux integrated over 
depth is given by equation (>1.31). In the relativistic limit we let fl -+ 1 and /?' — ► 1 
which gives 


Fr(E,n) 


f°° F 0 (E')dE' 
J b 


/jg(M) 

(1 -» 2 )e jb 

K»G(J*) 

(1 - ^E^iS - 1 ) ’ 


(34) 


where the second equality is for the power law injected flux. [Note that this expression 
may also be obtained by integration of equation (29) over r,.] For a given pitch angle we 
have a spectral index of (6 + 1) for this case. 


5) Total energy spectrum. Integration of F r over pitch angle (or F^ over depth) will 
give the total energy spectrum. However, if lim M — i G(fi ) ^ 0 (e.g. isotropic injection), 
the resulting expression diverges. This is because electrons with zero pitch angle never 
lose energy or change pitch angle; thus with a continuous injection there will be an 
infinite number of them from 0 < s < oo. This divergence disappears if as fi — ► 1, 
G(fx) -* (l-/i 2 )*, e > 0 and the total spectrum will be the same as that in equation (34). 
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The divergence will also be absent in the more realistic case of finite injection time or 
when collisions are included. 


C. Synchrotron and Collisional losses 


We need to consider both synchrotron and collisional losses when the ratio of these 
losses (iZ in equation (4)] is near unity. For non-relativistic electrons R ss 1 only when 
B is large. For relativistic electrons, however, synchrotron losses can be important for 
moderate values of B 2 if the density is low. There is no analytic solution for the general 
case including the synchrotron and collisional energy losses. Analytic solutions are 
possible for relativistic electrons because, as we have seen in sections A and B, Coulomb 
collisions and synchrotron radiation do not alter the pitch angles of relativistic electrons. 


1 ) Relativistic limit. In this limit we can ignore the last two terms in equation (1). 
The solution of this equation for uniform field and constant plasma density (i.e., E 1, 
dlnB/ds = 0 and constant R 0 ) is given by equation (A. 42) which reduces to 


*•(£,„, r e ) = F 0 {E.(E, r t )lG( M ) (1 + ( f + 


(35) 


where 


E*(E,n,r c ) = E 


1 + (£c/E) tan(r c //i£ e ) 

1 - (E/e c ) tan(r c / fi£ c ) J 


(36) 


and e c 2 = R 0 (1 — ^ 2 )- 

Note that equation (36) is valid only for E tan(r c /ji£ c ) < £ c . At a given pitch angle, 
F —> 0 at a depth given by 


< cr (E,v) 


Eljcr 

£c 


tan 1 



(37) 


where r, cr is defined in equation (30). As in the synchrotron-dominated case, for a given 
energy r' acr increases as fi increases, becoming infinite at n — 1. Then particles with high 
pitch angles are stripped away and the distribution narrows as depth increases. In the 
limit e c E (for R 0 1) synchrotron losses dominate and the flux and the critical 
depth reduce to equations (29) and (30) respectively. In the opposite limit ( R 0 < 1, 
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e c E) collisional losses dominate and equation (35) reduces to equation (14) as it 
must. 


2) Spatially integrated flux. We cannot integrate the flux given in equation (36) over 
pitch angle due to the complex /i dependence in E+ and e c but it is straightforward to 
integrate the flux over depth. We find 


= J™dsF(E,K,*) = (1 ^gr^f) /” Fo{E')iE' 

~ (1 + E 1 /«!)(<- l)E‘-‘ ’ 


(38) 


where the last relation is for the power law injected flux. In the two limits R 0 1 and 
R 0 < 1 this equation reduces to the expressions in equations (34) and (26) respectively. 


S) Total energy spectrum. If the injected distribution is narrow (i.e., <C 1) we can 
integrate equation (38) over pitch angle and obtain F TOT . We find 


Ke l/R.E*d 


F tot(E) - (£_!) a 2£6+l 


ln(H„.E 2 aJ) - 0.577 + f) ( ^ 

Jt=l 


2^ — fc 


(39) 


which reduces to 


TOT 


(E) 


K 


I (S-l)E 6 
< K 


— (1 - RoE 2 a 2 0 ) , 


R 0 E 2 a 2 o < 1 


ln(i2oJ5 2 a*/1.78) , R 0 E 2 a 2 0 > 1 


(40) 


{ (6 - l)R o a 2 0 E s + l 

Thus we have the expected spectral index for the collision dominated case at low 
energies and the index for the synchrotron dominated case (slightly modified by due to 
the logarithmic term) for high energies, provided that R 0 ot 2 ~ 1. This modification is 
due to the fact that collisional losses dominate for electrons with very small pitch angles, 
a 2 < 1 /R 0 E 2 . 


D. Non-uniform field (d\nB/ds ^ Q) 

Next we consider a non-uniform field for which d\sxB/ds ^ 0. We have no solution 
including collisional and/or synchrotron effects and a non-uniform field. [Ho (1986) has 
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given numerical solutions of the equations of motion for the case including synchrotron 
losses and converging fields, but he has not solved the kinetic equation.] A solution for 
the case with C = 5 = 0 in equation (1) was given in LP for the flux per unit area of 
the loop. Our solution, which is integrated over the cross-sectional area of the loop is 
different by a factor of B 0 /B(s), proportional to the change of the area with depth. In 
the absence of other effects B/(l — p 2 ) is a constant which leads to the solution 


F(E, 3) = F.(E)G (v/l-fl-^JB./BW) ~S^s) 


(41) 


which for G(p) = 2a~ 2 e ~ ( ' 1 becomes 


F(E,n,s) = 


2F 0 {E)B 0 

alB(s) 


exp 


(42) 


At any point s, the distribution has a dispersion given by a 2 0 B{s)jB 0 \ it is broadened 
by a factor of B{s)/B 0 . The flux at a = 0 is simultaneously decreased by a factor of 
B 0 /B(s), so that the number of electrons at a given depth integrated over pitch angle 
and area is constant with depth, as should be the case for zero energy losses. This result 
is independent of energy and therefore the demonstration by LP of the accuracy of the 
numerical code remains valid. 


1 ) Integrated Fluxes . We can integrate the solution given in equation (42) over pitch 
angle to obtain 

F,(E, s) = F,(E) . (43) 

It is clear that the energy dependence of the total flux F TOT will be the same as that for 

F a . 


III. Numerical Results 

We now describe results from numerical solutions of equation (1). To solve the 
equation we must specify the parameters of the background plasma (density and 
magnetic field) and the distribution of the injected electrons. We assume that electrons 
with the distribution given by equation (5) axe injected at the top, s = 0, of a symmetric 
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magnetic flux tube and solve the equation only for s > 0. We shall use both beam 
(<*x = 0) and pancake (<*j = tt / 2 ) distributions, but we note that the latter distributions 
may be inherently unstable and require acceleration perpendicular to the field lines. 
Electrons with fi < 0 at s = 0 are reflected back into the flux tube to simulate the 
symmetric geometry. Thus the total flux at s = 0 is equal to F 0 (E, /j.) + F(E, — /x,0). The 
knowledge of the geometry of the flux tube is not necessary here but will be essential for 
the evaluation of the angular and spatial dependence of the emitted radiation. 

For our purpose here all we need are the values of the coefficients d\n B / ds, S and 
C, which can be obtained from the variation of the density and magnetic field with 
depth s, B(s) and n(s). If B and n are constant only five constants, 6 , ai, a. 2 0 , n 0 , and 
R 0 = S/C , are needed for the solution. Unless otherwise specified we assume a fully 
ionized hydrogen plasma with £ = 1. 

A. Uniform Density and Field 

The parameter which determines the behavior of the electrons here is the quantity 
R defined in equation (4). This ratio depends on R 0 and the energy and pitch angle of 
the electrons. We describe the effects of all these parameters by considering plasmas 
with different values of R 0 and discuss the spectral and angular distributions at different 
H and E, respectively. We shall limit our discussion to electrons with energies between 
10 keV and 100 MeV. 

1. Collision Dominated Models 

For R 0 < 1 so that RJ3 3 7 2 <C 1 even for the highest energies of interest 
collisions dominate. Numerical results for non-relativistic electrons were given in LP and 
Leach (1984) and will not be reproduced here. We simply point out that the dispersion 
in pitch angle depends on the ratio r c /t] « t c /E 2 for energies < 300 keV. 

For higher energies such a simple scaling is no longer valid but a good 
approximation to the flux is given by equations (15) and (16). This is shown in Figure 2 
at r c = r) for a model with R 0 = 3 x 10 -9 , <5 = 5 and al = 0.40 for 16 keV, 300 keV, 

1 MeV and 10.6 MeV electrons. The solid lines are the numerical results for the pitch 
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angle distribution of the electron flux F(a)/F(a = 0) and the o’s denote the analytic 
results (eqs. [15] and [16]), i.e., for the same model in the limit R 0 = 0. The pitch angle 
distributions broaden with increasing depth due to diffusion except for E = 10 MeV 
where diffusion is small and almost overshadowed by the higher order pitch angle term 
discussed in Section II. 1. For all energies, the distributions are close to the analytic 
approximations even for angles near tt/2. 

In Figure 3 we show the same curves for the injected pancake distribution oi = 7 t/ 2 
and a 2 = 0.40. Here we have no analytic approximations. As shown by this figure, 
the pancake character (i.e., the maximum at ir/2) is lost very quickly for low energies 
(at r c = rj). Even for the highest energy shown (10 MeV) the maximum is shifted to 
a « tt/ 4 at r c = r). 

i) Spatially integrated flux, F r . In Figure 4 we show F r normalized to unity at zero 
pitch angle for energies 300 keV and 10.6 MeV. The solid lines depict beam injection 
(aj = 0) with a 2 0 = 0.40, and the dashed lines o x = ir/2 and = 0.40. Comparison 
of these curves with those in Figures 2 and 3 shows that for the beam injection, the flux 
at r c = rj is a good representation of the total angular distribution F r . For the pancake 
injection model the maxima at both energies occur at slightly larger pitch angles for the 
integrated fluxes than for the flux at r c = rj. 

ii) Pitch angle integrated flux, F M . In Figure 5 we plot F M versus normalized column 
depth t c /tj for 300 keV, 1 MeV and 10.6 MeV electrons for the beam injection model 
with a\ = 0.40 and <5 = 5. The solid lines are numerical results, and the o’s are the 
analytic results from equation (17). The agreement is good, and the behavior of F M is 
similar for all energies; for small values of t c /t), the flux is constant. As r c /r) increases 
beyond ~ 0.2 the flux falls off, and decreases as t~ s for r c /rj » 1. This behavior seems 
to be very general and fairly independent of the model parameters. 

2. Synchrotron Dominated Models 

As R increases, the synchrotron losses become more and more important and for 
large enough fields, they dominate even for non-relativistic electrons. As an example we 
eaxmine models with R 0 = 3 x 10 3 (or B 2 /n « 1), which corresponds to B = 10 4 G 
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and n 10 8 cm -3 , an extreme but unlikely condition for solar flaxes, or to B ~ 10 12 G 
and n ~ 10 24 cm -3 , which may be representative of conditions in the magnetosphere 
of a neutron star. Except for very small pitch angles, synchrotron losses dominate for 
the entire range of energies considered here ( R ~ 1 at 1 keV and a = 7t/4) and the 
distributions behave as we expect from the solution including only synchrotron losses, 
i.e., equation (A. 24), which reduces to equation (29) for E > 1. In Figure 6, we compare 
the numerical results (lines) with the analytic ones (o’s) for 300 keV and 10.6 MeV 
electrons at a depth of r, = 1 for a source distribution shown by the dashed line. The 
behavior of the non-relativistic electrons is qualitatively similar to that of the relativistic 
electrons; as the depth increases the distribution narrows and for a given depth, F — * 0 
at some critical angle ca cr = cos - 3 (/i cr ). [See equation (31).] Note that the numerical 
results do not quite agree with the analytic results on the value of the critical angle. 

This is due to the finite size of the angular grid used for the numerical calculations; 
the difference between the values of the critical angle is always one grid spacing and the 
numerical results axe inaccurate at the point next to the critical angle. 

i.) Very small pitch angles: At large depths most particles acquire very small pitch 
angles (a — ► 0) and the numerical results begin to deviate from the analytic solution 
beacuse the ratio R becomes less than one and collisions become important. This effect 
is shown in Figure 7 in which we plot the zero pitch angle flux (actually FE 6 ) versus 
energy for depths r, = 5 and r 4 = 10. The analytic results are shown by the dashed 
lines, and as we can see, for non-relativistic energies, the flux at zero pitch angle diverges 
as E — > 0. This can be deduced from equation (A. 24) which in the zero pitch angle limit 
reduces to 

F{E, r„ 0) = F 0 (E)G(fi)e 2r ^^ . (44) 

The angular width of the distribution decreases exponentially; 

tt ]«rV J e- 2r '^, (45) 

so that the flux integrated over pitch angles, oc F(a = 0)a 2 , is finite. 

The analytic solution breaks down not only because the effective value of the ratio 
R decreases as the critical angle approaches zero, but also because the gradient of the 
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pitch angle distribution increases, and collisions! diffusion can no longer be ignored. 

From equation (1) the ratio r D , of the diffusion term to the energy loss terms (both 
collisional and synchrotron), is of order 

r 0 = £E[/JVc 2 (1 + f^R.a 2 )]-' . (46) 


so that for pitch angles less than 


Or, = 


\/4/3JE + 1 — 1 


|!/2 


(47) 


r D exceeds unity and the analytic results are not valid. The depth at which a 2 D 
thus 



4 

0 2 7 2 <*l 


= a 2 r is 
(48) 


For depths beyond t 3 j we do not expect the analytic solution to be useful. Table 2 gives 
values for a 2 and t s< j for various energies and R 0 = 3000. 

Note that a 2 D increases and t 3C i decreases with decreasing energy. For 16 keV 
electrons, the analytic solution is incorrect even at r,= 1, while for 10 MeV electrons 
the solution may be used for large values of r,. This behavior is evident in Figure 7. At 
t s = 5 and 10 the numerical solutions fall away from the analytic solutions for energies 
less than 2 MeV and 5 MeV respectively. 


S. Models with Intermediate R 0 

Next we consider a model with R 0 = 1.3 ( B 2 /n ~ 4 x 10" 4 ) so that synchrotron 
losses are important for energies greater than the rest mass energy, and collisions 
dominate for E < 300 keV. For low energies the scaling of the pitch angle distribution 
with tJE 2 described in item 1 above is still valid. For relativistic electrons with 
energies > 5 MeV the synchrotron energy losses dominate and we obtain distributions 
similar to those expected from equations (29) and (30). According to these equations, 
the pitch angle dispersion will be the same (or curves of F versus a will have the same 
shape) for depths r, oc 1 /E- We have found that the numerical results agree with the 
analytic values of (eq. [30]) to within 10 percent for r t E = 1. 


17 



In Figure 8 we plot the integrated flux distributions F r /F r (a = 0) for the beam 
(solid lines) and pancake (dashed lines) distributions with a 2 0 = 0.40 at 300 keV and 
10.6 MeV. The curves axe intermediate between those given in Figures 4 and 6 with the 
effect of synchrotron losses evident 10.6 MeV and the effects of collisions dominant at 
300 keV. Note that the pancake injection model still has a non-zero maximum, but very 
few electrons with a > ir/2. 

In Figure 9 we plot F <J (E,r), the flux integrated over pitch angles, versus energy 
for the beam and pancake distributions with aj = 0.40 at r, = 0.1 and 1.0. At 
larger depths the curves fall off at low energies due to increasing collisional losses with 
decreasing energy and they steepen at high energies due to increasing synchrotron losses 
with energy. The synchrotron losses for the pancake distribution are more noticeable 
because a larger proportion of the electrons have high pitch angles and therefore higher 
synchrotron losses. 

4- Spectral index and curvature 

We now discuss the spectral index of the spatially integrated flux F r (E, n). In 
Table 3 we give the spectral indices, slopes of power law fits of F r versus E curves, for 
different pitch angles in different energy ranges for R a = 3 x 10 -3 and R 0 = 1.3. For 
energies < 300 keV, we find essentially the same angles and models with spectra for 
R a = 1.3 in general steeper. Here collisions dominate and the slope is approximately 
(5 — 1. For relativistic energies (E > 10.6 MeV), we find different behavior for the two 
different cases. For the beam injection and a < tt/ 2, the slope is slightly larger than <5 — 1 
for small R, and nearly 6 + 1 for large R and nonzero pitch angle. These values are in 
good agreement with the results expected from the analytic calculations of Section 2.3. 
For a > 90°, however, this slope becomes enormous, reflecting the absence of collisional 
diffusion for high energies. For uniform fields, very few relativistic electrons are scattered 
to high pitch angles. 

For the pancake injection, the spectral indices are larger for relativistic energies. 

This is again due to the larger proportion of electrons with pitch angles near a = tt/2. 

To show the effect of synchrotron losses, in Table 4 we give spectral indices fit to 
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the total integrated flux F TOT (E) for <5 = 5, a 2 0 = 0.40 for both pancake and beam 
injection, and for R a = 3 x 10 -9 , 3 x 10 -3 , 1.3, and 3 x 10 3 . We see the expected 
results for the beam distributions; for low energies (16 to 300 keV) the index is ~ 6 — 2 
for low fields and increases for high fields. For high energies (10 to 76 MeV) the index is 
~ 6 — 1 for 5 = 0 and increases to slightly less than 6 + 1 for B > 2000 G. Thus our 
analytic results give accurate answers for the total flux, which is important for spatially 
unresolved sources. 

For the pancake injection, the results are the same for the non-relativistic and 
semi-relativistic ranges, but as was the case for the other integrated fluxes, the spectral 
indices are larger for relativistic energies. 

B. Uniform Density, Converging Fields 

We now turn to models which have non-uniform (in particular, converging) 
magnetic fields. This produces two changes. First, since dhiB/ds is not zero we must 
include the second term in equation (1). Second, the coefficient 5 varies along the field 
lines so that R 0 is no longer constant. We parametrize these models with B 0 , the value 
of the magnetic field at s = 0, and the parameter Sb = dlnB/ds which we assume to be 
a constant; b(s ) = B{s)/B 0 = exp(sSfc). The effects of convergence become important 
when Sb is of order of the coefficients C or 5 for collisional and synchrotron losses, or 
when the dimensionless column depth rj = J Sbds , which in this case is equal to In b, is 
greater than r c or r 4 . Since we have assumed constant values for and C the relative 
importance of these two effects remains constant along the field lines but varies with 
the energy and pitch angle of the electrons, convergence becoming more important for 
higher energies and larger pitch angles. The coefficient S , however, varies along the field 
lines (S = S 0 b 2 ) so that the relative importance of the synchrotron losses increases with 
depth and, of course, with energy. Thus at sufficiently large values of depth s or energy 
E, t , = {S 0 /Sb){b 2 — l)/2 exceeds Tj, and synchrotron losses become dominant. 

1. Strong convergence and loss-cones 

An important effect of the field convergence (when Sb > C and 5) is to produce 
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pitch angle distributions with distinctive “loss-cones” in the direction of the field lines. 
We have no analytic solution which accounts for convergence along with snychrotron 
radiation and collisions so we present numerical results only. 

To illustrate the interplay between these processes we consider a model with 
isotropically injected flux (i.e., = oo), with C = 2 x 10 -13 cm -1 , S 0 = 6.5 x 10 -16 cm -1 

(e.g., B 0 = 100 G, n 0 = 10 10 cm -3 ), and Sj, = 9 x 10~ lo cm _1 , which is large enough to 
insure that the effects of convergence dominate over synchrotron and collisional losses for 
150 keV < E < 100 MeV. The pitch angle distributions at the top of the loop for 16 keV, 
1 MeV and 10.6 MeV electrons axe shown in Figure 10 for a model with finite length 
■Smax = 2.4 x 10 4 km or r 6imax = 2.3. Electrons reaching s > s max are ignored. In 
the absence of collisional diffusion electrons with initial pitch angles o, n < a cr with 
sin a cr = fw£ 2 reach s raax and escape. The 10.6 MeV curve in Figure 10 shows this 
rapid decrease in electron flux a < 27° and a ir — 27° . 

This effect is less pronounced at lower energies due to collisional diffusion (as shown 
by the 16 keV curve in Figure 10) and also for large magnetic fields and high energies 
due to the effectes of synchrotron losses. This is shown by the dashed line in Figure 10 
for 10.6 MeV but with 5 increased by a factor of 40. 

2. Weak convergence, relativistic electrons 

In some situations the presence of convergence can have significant consequences 
even when convergence is not dominant. As an example we consider collision-dominated 
models with C = 2 x 10“ 6 cm _1 and 5ft = 3.7 x 10“ 9 cm _1 . In spite of the fact 
that 5* <C C the effect of the convergence term is approximately equal to that of the 
pitch angle diffusion term for relativistic electrons, and can have a significant influence 
in broadening the pitch angle distribution when collisions are ineffective in doing so. 

In Figure 11 we compare the distributions of the uniform (U) and converging (C) field 
models with isotropic injection, E = 10.6 MeV, r c = rj and two values of R 0 . 
Synchrotron losses dominate for R 0 « 1 so that there are few electrons with high 
pitch angles and the convergence has essentially no influence. For ^ < 1, however, 
there axe some electrons with high pitch angles at this depth and the effect of field 
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convergence is to increase the number of such particles even further. This effect is 
particularly noticeable at higher energies, but not too high an energy when synchrotron 
losses dominate. Even a small amount of convergence results in a substantial number of 
these reflected electrons as shown by the increasing divergence of the C and U curves in 
Figure 11 for pitch angles a > 7r/2. 

To further illustrate this effect, in Table 4 we give the values of spectral indices from 
power law fits to the spatially integrated flux F r for the cases with small convergence 
discussed above. The values of the spectral index m x for low energies 16 to 300 keV 
are similar to those of the uniform field models. At higher energies, convergence has 
little effect for a < tt/ 2, and the spectral indices are again similar. But for a > 7r/2 
the indices are much smaller for the converging field cases since there axe more reflected 
electrons. 

C. Solar Flare Models 

It is believed that during the impulsive phase of a solar flare high energy 
particles (10 keV to > 10 MeV) are accelerated in a coronal magnetic loop of length 
10 9 to 10 10 cm, with density n as 10 10 cm -3 and a magnetic field of a few hundred 
gauss (Kundu 1983, Lu and Petrosian 1989). Below the transition region the density 
increases rapidly. It is not known how the magnetic field vanes with depth but it is 
suspected to increase to a few thousand gauss at the photosphere. Above the transition 
region the plasma is fully ionized (C = C' or f = 1) but below the transition region the 
temperature decreases, and the atmosphere becomes neutral. 

For n = 10 10 cm -3 and B 0 = 100 G, the quantity R a = 3 x 10 -3 at s = 0 and 
could be as large as R 0 = 0.3 at the transition region with r c = 0.0014 and B tI = 10 3 G 
(or Tfc(s tr ) = 0.07). Thus convergence may be important but collisions! and synchrotron 
losses are negligible except for E < 20 keV and E > 50 MeV, respectively. Below the 
transition region density increases rapidly and unless the magnetic field scale height 
begins to decrease rapidly synchrotron and field convergence effects become negligible. 

In Figure 12 we plot the pitch angle distributions for 300 keV and 10.6 MeV electrons at 
the top of the loop (solid lines), and at r c = r) (dashed lines) for models with isotropic 
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injection, <5 = 5 and R 0 = 3 x 10 -3 ( B 0 = 100 G and n = 10 lo cm“ 3 ) at s = 0 which 
increases to R 0 = 0.3 at the transition region. As mentioned above, the convergence will 
be important in the corona and at the top, we have loss-cone distributions. Below the 
transition region collisions dominate and we see typical behavior at both energies. 

Not surprisingly, the integrated fluxes behave as in the uniform density collision 
dominated cases. Unless the magnetic field is very high above the transition region, 
there axe no discernable effects of synchrotron losses on the integrated fluxes. Since 
T c (s u )/rj <C 1 for most of the energy range, the integrated fluxes F r and F TOT are 
dominated by the behavior in and below the transition region and the integrated fluxes 
look like those for collision dominated models. 

One quantity which is readily available from observations is the energy spectral 
index. To show the effects of various parameters on this quantity, in Table 5 we give 
values for spectral indices m 7 (for energies between 10.6 and 76 MeV) and m x (between 
16 to 300 keV) of the integrated flux F r for models with 6 = 5, including: uniform field 
models with beam and pancake injection with a 2 = 0.40 and R 0 = 3 x 10~ 3 and 1.3, 
models with isotropic injection, B 0 = 100 G ( R 0 = 3 x 10 -3 at s = 0) and mirror ratios 
of b tI = 1, 2, 5, and 10; strongly beamed models with a 2 = 0.04, B a = 100 G, and 
b« = 2 (to show how the effects of small convergence combine with narrow beaming); 
and models with a 2 — oo, B 0 — 2000 G (.Ro(top) = 1-3), and bt T = 1 and 10. For all the 
weak field cases, even for low levels of convergence, and for the sharply beamed injection 
the values of Am = m 1 — m x are close to zero for a > n/2. For the high field case 
the value of Am is larger for a > tt/ 2 and nearly as large sis for the uniform weak field 
case. Thus we find that a high field (here B(s u ) = 20,000 G) can negate the effects of 
convergence. 

Values of the spectral indices for the total integrated flux for these models are given 
in Table 6. For the cases with B 0 = 100 G we find that F ror (E) behaves as in the 
collision dominated uniform field cases, since F VOT is dominated by the behavior of 
the distribution below the transition region, where the collisional scale length is much 
smaller than the scale of convergence. For the model with B 0 = 2000 G increased 
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synchrotron energy losses lead to larger slopes in all regimes. 

D. Gamma-Ray Burst Models 

Gamma- Ray bursts are believed to occur in the magnetospheres of neutron stars 
with magnetic fields ranging from 10 12 G (as deduced from absorption lines) to less than 
10 11 G (Matz ei al. 1985) to allow high-energy (E > few MeV) 7- rays to escape without 
pair production. The power law energy spectra observed in some bursts to energies 
greater than 30 MeV are indications of the presence of non-thermal electrons accelerated 
to similarly high energies. These conditions require the examination of the transport of 
the relativistic electrons in variable magnetic fields and possibly inhomogeneous plasmas 
reaching densities of as high as lO 3 *’* electrons/cm 2 at the surface of the neutron stars. 
All three effects on the electron transport of the relativistic electrons considered here 
(synchrotron, collisions and field convergence) may play significant roles in this process. 

If the electrons are accelerated in the magnetosphere, where R 0 is expected to be 
very large, they suffer synchrotron losses. They quickly lose all of their perpendicular 
momentum and slide along the field lines with near zero pitch angles. In this phase 
the synchrotron loss formulae provide an accurate description of the problem up to 
some column depth. As shown in Section III, for large values of t s collisional diffusion 
becomes important. Since the field lines are anchored to the neutron star, the electrons 
hit the surface where they undergo pitch angle diffusion and emit more synchrotron 
or bremsstrahlung radiation if the density is very high. Equations dealing with both 
synchrotron and collisional processes will be applicable. It is unlikely that field 
convergence will play a significant role dumig this transport because the electrons will 
have small pitch angles and the scale height of the convergence S 6 -1 (of the order of the 
neutron star radius) is probably much larger than the scales C~ l or S~ l for collisional 
and synchrotron process. 

One can consider the problem in two steps, starting with the magnetospheric part. 
Consider a loop structure in the magnetosphere of a neutron star, with a magnetic field 
of 10 11 G and an ambient electron density of < 10 24 cm“ 3 . For this case R 0 > 32.3, 
and synchrotron losses are dominant. The scale length for synchrotron processes is 
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S -1 = 1.55 x 10 -3 cm, which even for a length of about one cm gives t s > 10 3 which is 
much larger that T„d (eq. [48]) where diffusion must be taken into account. [See Table 3.] 
We cannot calculate accurate analytic or numerical solutions for this case. The analytic 
solution diverges and the width becomes very small. We may, however, obtain a good 
approximation to F^, the flux integrated over pitch angle, using the analytic solution 
with synchrotron losses only. Collisional diffusion has no effect on the value of F t 
therefore we expect this calculation to be useful even for extremely large values of r s . 
This flux asr, -too can then be used as input for the second part of the problem, where 
we solve the full equation starting at the surface of the star. 

For large t s the electron distribution is narrow and we may use the small angle 
approximation to the synchrotron-only solution. We have from equation (A. 27), 

X = 1 - 2a 2 /a 2 cr , (49) 

where a cr is given by equation (>1.29). For isotropic injection, the solution 
equation (A. 24) for the electron flux for large t„ becomes 

F(E, a, r,) = F 0 (E)e 2r ‘' l} ‘ f e- ai ' a * , (50) 


where the dispersion is given by 


2[(£-l)7 2 +«7 + 3] 


a. = 


(51) 


The width decreases and the flux for zero pitch angle increases with t b but the flux 
integrated over pitch angle is finite and is given by 


F lt (E,T,-+oc) = F 0 (E)[(6-l) 7 2 +6y + 3]- 1 . (52) 


If we include collisions the flux at a = 0 does not diverge, and the width remains 
finite but small. Since the width is too small to include on the pitch angle grid which we 
use for numerical solution, for the calculation of the flux below the surface of the star we 
inject all of the flux given by equation (52) at a = 0. 
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The plasma at the surface is not an ordinary plasma. The electrons may be 
degenerate and there are heavy ions. Since the ratio £ of the diffusion coefficient to 
the energy loss coefficient is highes by a factor equal to the atomic number of the ions, 
we have the situation in which synchrotron energy losses and collisions! diffusion are 
important. The pitch angle distribution of 1 MeV electrons for a model with a density of 
Fe ions of 10 27 cm -3 and a magnetic field of 10 12 G (or R a = 1.2) is shown in Figure 13. 
Each curve is marked by the value of r c ~ t, for that depth. The flux at r c = 0 is a 
delta function, but the distribution is nearly isotropic even for small depths. Since this 
is true, when looking at any emission the important quantity is the total electron flux 
F tot (E). The spectral index at the top of the magnetospheric loop is 8 = 3 and the 
spectral indices for the delta function flux at the surface are m x — 2.3, m m = 3.4, and 
m 7 = 4.7. The indices for the total flux are m z = 2.6, m m = 4.4, and m 7 = 5.9. 

As we can see, this combination of effects (collisional diffusion plus synchrotron losses) 
results in steepening of the spectrum at every energy, and large breaks in the spectra 
between non-relativstic and relativistic energies. Such breaks have been observed in the 
spectra of 7- ray bursts. A detailed comparison with observations is beyond the scope of 
this paper and will be treated in a subsequent work. 

IV. Summary 

We have extended the Fokker-Planck analysis used in LP to include ultrarelativistic 
electrons and effects due to both collisions and synchrotron emission. We have solved 
the resulting kinetic equation (eq. [1]) analytically in some simple limiting cases, and 
have shown that the numerical results for these cases agree well with the analytic 
results, thus verifying the accuracy of our numerical code. We have shown that 
the analytic results provide a useful guide for the qualitative description of the flux 
distribution of the non-thermal electrons and examined the effects of non-uniformities 
in the magnetic field, in particular converging magnetic field geometry. 

Some of the features of our results axe: 

(1) In situations dominated by collisions equations (15) and (16) give a fairly 
accurate description of energy spectra of both relativistic and non-relativistic electrons 
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and of the flux distribution for the gaussian pitch angle distribution of accelerated 
electrons. Even though these equations axe small pitch angle approximations the are 
good for much larger angles than expected. Equations (17), (24) and (28) give the flux 
integrated over pitch angle, depth and both, respectively. 

(2) For problems where the synchrotron process is dominant we have obtained a 
complete time-dependent solution (eq. [>1-21]) and a steady state solution (eq. [A.24]). 
The relativistic limit of the latter acquires a simple form given by equations (29), (32) 
and (34). 

(3) If both synchrotron and collisional losses are important then Equations (35), 

(38) and (39) give the fluxes in the relativistic limit. Analytic solutions are not possible 
for the less likely situation of non-relativistic energies with both synchrotron and 
collisional losses. 

(4) The above results Eire for uniform magnetic fields or are applicable when 
the magnetic field variation scale is much larger than the scales for collisional and 
synchrotron losses. For cases when this is not true one must resort to numerical 
solutions as described in the text. 

The discussion of the behavior of non-thermal electron distributions included in this 
paper will be important for our future work, in which we will use the numerical solutions 
to calculate the expected radiation during the impulsive phase for solar flares or cosmic 
gamma-ray bursts. Examples describing models for these bursts will be presented 
and we shall compare the radiation signature of such models with observations and 
determine the distribution of the accelerated electrons responsible for those radiations. 
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Appendix A: Analytic Solutions 

We now present some approximate analytic solutions of equation (1). 


1.) Collision Dominated Solutions 

Consider equation (1) with d\nB/ds = 0 and 5 = 0. We have (with dr c = Cds) 


^d7 c ~ dE 


±± ' 
E 2 dp . 


n 


(A.l) 


where £ = C' /C(= 1 for a fully ionized background plasma). The small pitch angle 
solution for non-relativistic electrons is given in the text. This solution (eq. [8]) cannot 
be extended to high energies because as E increases the term ignored (of order a 2 d$ds) 
becomes more important relative to the diffusion term. For extreme relativistic electrons 
we may ignore the pitch angle diffusion term on the right hand side of equation (A.l), 
and we have 


^!h c ~dE~°' 


(A.2) 


which for injected $(£,/*, 0) = $ 0 (E)G(n) has the solution 


$(E, n, t c ) = $o{E + r e /fi)G(ti) . 


(A3) 


For relativistic energies, fi = 1 and F = c$. In the small pitch angle regime and for a 
gaussian pitch angle distribution G(fi) = 2a~ 2 e a / a °, equation (A.3) gives 


F(E, n, r c ) = F 0 (E + r c ) 


2 e -°V^ r 6 e r c a 2 

a 2 [ 1 + 2(£ + r c ). 




(A.4) 


where 6 e = d\nF 0 (E)/d\nE\E+ rc and is equal to the spectral index when F 0 is a power 
law. Note that the term in the brackets which is of order 6 e a 2 has an effect which is 
opposite of the diffusion term; namely that dispersion decreases with increasing depth 
r c . 

For lower energies or smaller at 2 the diffusion term becomes important and in the 
intermediate energies it increases the dispersion with depth at a rate determined by the 
value of l/Ea 2 . To account for this effect we treat the pitch angle diffusion term as a 
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perturbation. The fractional perturbation we denote by <f>\. In the small pitch angle 
regime <j>\ obeys the equation 


where 


the solution of this is 


d<p\ _ d(f>i 4£ / a 2 - «r \ 

dr c dE ~ E 2 ot\ \ aj J ' 


2 _ a o 

’ 1 + Sal/2(E + T c ) ’ 


, _ 4 if a 2 -a 2 z \ Tc 

<* 2 * jE(E + r t 


) 


(A.5) 


(A.6) 


(A.7) 


Combining this with equation (A. 4) we then obtain, in the small pitch angle regime, 


F{E, n, t c ) = F 0 (E + r c )- 


ai 


where 


and 


a] = a] + [*r JE(E + r«)] , 




a; 


(A.8) 


(A.9) 


(A.IO) 


[1 + 6 e r c al/2(E + r c )] ' 

This form is similar to the non-relativistic solution (equations [8] through [10] in the 
text). With a small modification of a 2 to the form 

.2 


Ck 

a<: ~ [1 + S e r c a 2 c / 2(2 + rj + r c )] ’ 


(All) 


with a 2 given by equations (9) and (10), and addition of the velocity factors as in 
equation (8) we can combine the solutions of the two limiting cases as 


m y 

mv + r c))J 

( m y 

\mn + Tc))J 


F(E,h,t c ) 


■ 


— — — -F 0 [E(ri + r c )][l + S e T c a 2 / 2(2 + rj + r c )] 

Q : c 

2 e -a i /&l 

2 F 0 [E(t, + t c )} , (A. 12) 

a c 


which describes the solution with a high degree of accuracy and to much larger values of 
the pitch angle than expected considering the small pitch angle approximation. 
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2.) The Uniform Field Synchrotron Solution 


It is well known (see e.g., Jackson 1962) that the rate of change in /?u = the 
velocity parallel to the magnetic field, due to synchrotron losses is zero; i.e., (/3||), = 0. 

If collision al losses are neglegible and the B field is uniform, it is then convenient to 
rewrite the kinetic equation in terms of and /? x = /?(1 — /x 2 ) 1 ^ 2 . Furthermore, since /3\\ 
is constant the time dependence can easily be incorporated and equation (Al) becomes 

• ( - 413) 

The function / is now the distribution function in terms of and we can use the 

Jacobian of the transformaton from to (/?,,/? l) to give us the relation between / 

and /, the distribution in terms of ( E , ^); 


[i 3 (i-j< ! ) I/: 7(E,A)] 


(-414) 


The tilde denotes quantities which axe to be functions of (/3 |,/?l); i.e., 7 = (1 — /? 2 — 
E = 7 — 1, and fi — + /^l) 1 ^ 2 - If we make the substitution 


h = /?x(l - A ~ % ) 1/2 f = 7(f - 1) 1/2 (1 - A) , (44-15) 


and define y = ctS and r, = sS equation (A. 13) becomes 


dr t 

We define a new variable u by 

du = 




401 


0.(1 - 0 l ~ 0?) ,/2 ’ 


so that 

u = — 7„sech -1 (7|^x) and = — sech(-u/7 ( ) , 

7| 

where 7 = (1 - /? 2 ) -1/2 . Equation (A. 13) is simplified to 

dh dh _ dh 
d^ + ^dr t ~ du ' 


(A.16) 


(A.17) 


( 44 . 18 ) 


(.4.19) 
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If at r, = 0 the distribution is given by g(y)f 0 (E, n) so that 


&.(«,0,y) = [i(7 2 - 1) 1/2 (1 - ?)ME,ii) ] 

. L -I E(u),n(u) 


(A. 20) 


Equation (A. 13) has the simple solution h(u,r,,y) = g(y — T 3 /0^)h o (u + r a //?u), and the 
distribution function in terms of energy and pitch angle / = h(u,r 3 ,y )/ - /z 2 ) is 


Wm.») = [*? - D' /2 (i - m^)] £( . + ,. 




(A.21) 


It is easy to show that 


E(u + r ,/$,) = £(X) = (7, - X)/X 

«» + = A(X) = /J,T P /(T, 2 - X 2 )'/ 2 , (A. 22) 


where 




/Al7||) 

/At.) 


(A. 23) 


Returning to the time independent case of continuous injection, for which g(y) is a 
constant, we find for the flux F = (icf 

7 ?F(£(X) 


F(E,h,t 3 ) = 


7 2 (1 


:),/i(X),o) / i - x 2 
-/i 2 ) \X 2 (7 2 — X 2 ) 


(A.24) 


where F(E,n, 0) is the injected flux. 

At t s = 0, X = 7 ( /7, with increasing r,, X and the flux F decrease monotonically 
and become zero at a depth given by 


T»cr = 



(A. 25) 


The value of T scr decreases with increasing energy; particles with higher energy travel 
shorter distances, and for a given energy t 3CT is smallest at n = 0; particles with high 
pitch angles lose energy more quickly than those with smaller pitch angles. 
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At a given depth r 4 , the flux approaches zero at a critical value of the pitch angle 
a cr = cos -1 (^i cr ) which we may calulate by solving the trancendental equation 


tanh 


-l 


~ r ) 1/2 
fl/Ac r 


(A. 26) 


M l-^ 2 ^r) 1/2 J 

For large depths r, — ► oo, /i cr — » 1, and the distribution becomes infinitely narrow, and 
at the same time, the flux at fi = 1 gets infinitely large. For p — ► 1, or a 2 — ► 0, we have 

02^2 a 2 e 2r, / Py 


X = 1 - 


(A. 27) 


E(X) = E, £(X) = n, and 


F(E,n = l,r a ) = F(E,l,0)e 2r -'^ 


(A. 28) 


For a cr < 1, we may solve equation (>1.26) and find value for the critical angle given by 


a 


2 _ 
cr 


4 e 2 r./^ 7 

0 2 7 2 


(A. 29) 


Even though the flux at a = 0 diverges for large depths, the width becomes zero, so that 
the flux integrated over pitch angle, which is of order F(fi = l)a 2 r , remains finite. 

The solution has a particularly simple form for ultrarelativistic electrons. We let 
(3 — ► 1 in equation (A.21) and find X = 7|(1 — T,/r, cr )/7 ■C 1, ^i(X) = and if 
F(E,n,0) = F 0 (E)G(n) then 


F(E,n,s) 


Fq[E/( 1 ~ -r,/T Jcr )]g(/x) 
(1 - r,/r, cr ) 2 


(A.30) 


where r 4cr = p/E( 1 — /i 2 ). This solution is discussed in Section II. 

In Section II we discuss the flux integrated over depth for given pitch angles. 

This would pertain to the case for which the field lines are straight and for spatially 
unresolved observations. We may obtain the flux integrated over depth either by 
integration of the solution in equation (A.30) over r, or by solution of the full equation 
(eq. [A. 13]) integrated over r s . We find 


F r (E,n) 



,H,T,)dT a 


3u 

7 2 (1-^ 2 ) 



dE-^la 



(-4.31) 
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where we must keep /3 ( | constant in the integrand. The integration is complex for the 
general case due to the complex dependence on /?' in the integrand. For an initially 
isotropic distribution, G(0p/0') is constant and F s oc j3'~ 2 F 0 (E')dE ' . 


3.) The Relativistic Energy Loss Solution 

As in the collisional case, the synchrotron pitch angle change term is much less 
than the energy loss term for E > 1. The relativistic electrons lose energy at constant 
pitch angle until they become non-relativistic when the collisional diffusion and p, terms 
become important. Here we ignore these terms, keeping in mind that we can add those 
terms as perturbations as in item ( 1 .) above if neccesary. We take the limit j3 — *■ 1 and 
let d\n.B/ds = 0 in equation (1). If we define e 2 = CfS{ 1 — /i 2 ), then 

/*!£ = (1 + E 2 /a|§ + . (A.Z2) 

and if dt] = dE/( 1 + E 2 /e 2 ), then 


rj - e c tan -1 (£/£<;) , E(tj) = e c tan(r?/£ c ) , (A.33) 


and 


5$ d$ 2 A 

= - 3 — H tan 

5r c drj e c 


(£>*• 

Next we let $(rj,n,T c ) = <£(/>)t/>(r?) where p = r) + r c //i so that 


dr c 


d<j> d* 
„ ” ^ dp ’ dr} 


1 ^ 


and equation (A. 34) reduces to 


drj> 2 

— tan 

dr} £ c 


fe) 


0 


which has the solution 


This gives 


V>(»7) = ^ 2 ( t // £ c) 


$(r?,/i,r c ) = <f>(p)ip 0 sec 2 {r}/e c ) , 


(A.34) 
(A. 35) 

(A. 36) 

(A. 37) 
(A. 38) 
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where rp 0 is a constant of integration. 

For the initial condition 

$(£, p, 0) = 4>{p = r))\l>{r}) = c " 1 F 0 [E(rj)]G(p) 


we then find 


<t>(ri)tl>o — c 1 G{p)F 0 [e c tan(r// e c )] sec 2 ^/ e c ) . 


Thus 


and, finally 


4>(p) = 


F(E , r c ) = c$(F, /x, r c ) = 


F 0 [g c tan(^>/ e c )]G(/i) sec 2 (p/e c ) 

crp 0 


c<t>(p)ipo Fp [e c tan(/?/e c )]G(/i) sec 2 Q/ g c ) 


sec 2 (ri/e c ) sec 2 (r}/s c ) 

We substitute for p and 77 to obtain the flux int terms of E and p; 

F(E,u,n) = F.{£.(E, x, r«)]C ( > .) (1 + ( f j 


where 


E m (E,p,r c ) 



+ g c tan (r c /^e c ) 
- Ftan(r c /^e c ) 


(A.39) 

(A.40) 

(A.41) 

(A.42) 

(A.43) 
(A. 44) 
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Table 1 


A comparison of analytic widths with numerical results for 10.6 MeV electrons, 
with 8 = 5, = 0.04, 0.10, 0.40, and 40, at depth r c = rj. The subscript * denotes a 

numerical result. 


4 /Ea 2 0 6aH 2 

* 2 0 

<* 2 c(v) 

a 2 c(v) 

al(ri) 

5.00 0.10 

0.04 

0.13 

0.12 

0.11 

2.00 0.25 

0.10 

0.20 

0.16 

0.17 

0.50 1.00 

0.40 

0.50 

0.31 

0.31 

0.005 100.0 

40.0 

40.10 

0.81 

0.61 


Table 2 



Critical values for the dispersion a 2 p and depth 



E 16 keV 

300 keV 

1 MeV 


10.6 MeV 

a\ 0.09 

0.01 

3.50 x 10' 

-3 

1.80 x 10~ 4 

T »d 0.83 

3.40 

7.00 


42.0 
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Table 3 


Spectral indices of the spatially integrated flux F r for energy ranges 16 to 300 ke\ , 
600 keV to 5.32 MeV and 10.6 to 76 MeV, respectively. 


a 

Ro = 3 

m x 

x 10" 3 
m m 

m 7 

R 0 = 

m x 

1.3 

TTlm 

m 7 

Beam, 

a 2 = 0.40 






54° 

4.00 

4.20 

5.62 

4.16 

5.68 

5.99 

126° 

4.23 

6.86 

14.7 

5.17 

14.3 

15.9 

Pancake , a 2 = 0.40 






54° 

3.98 

3.99 

5.55 

4.14 

5.46 

5.99 

126° 

4.15 

6.41 

20.3 

4.88 

17.4 

19.6 

Isotropic, a 2 = oo 






54° 

3.98 

4.03 

5.57 

4.13 

5.50 

5.99 

126° 

4.18 

5.87 

24.7 

5.00 

17.6 

19.6 

Isotropic, ag = oo, 

convergent 




54° 

3.98 

3.99 

5.06 

4.03 

5.18 

6.00 

126° 

4.17 

5.72 

11.2 

4.39 

15.5 

14.6 
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Table 4 


Spectral indices of F TOT for 8 = 5 and a \ — 0.40. 


R 0 

Beam 

m x 


m 7 

Pancake 

m x 

m m 

m 7 

0 

3.20 

3.91 

4.05 

3.22 

3.90 

4.05 

3 x 10" 3 

3.20 

3.94 

4.92 

3.22 

3.97 

5.54 

1.3 

3.32 

4.86 

5.80 

3.40 

5.38 

5.97 

3 x 10 3 

3.50 

6.11 

6.02 

4.38 

5.83 

5.98 
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Table 5 


Spectral indices of F r for different pitch angles for various flare models. The 
quantities m x and m 7 vare as defined for Table 3. 



a = 

54° 

a = 

126° 


<*1 

B 0 (G) 

b 

m x 

m 7 

m x 

77T 7 

0.40 

0 

100 

1 

4.05 

4.06 

4.89 

10.30 

0.40 

7t/2 

100 

1 

3.96 

4.02 

4.54 

9.89 

0.40 

0 

2000 

1 

4.04 

4.20 

4.89 

12.70 

0.40 

• 7r/2 

2000 

1 

3.96 

4.18 

4.54 

12.8 

0.40 

0 

10 4 

1 

4.03 

9.62 

4.83 

24.00 

0.40 

7t/2 

10 4 

1 

3.98 

10.0 

4.53 

26.30 

oo 

— 

100 

1 

3.98 

4.02 

4.65 

9.93 

oo 

— 

100 

2 

3.94 

4.05 

4.39 

5.02 

oo 

— 

100 

5 

4.04 

4.06 

4.43 

5.00 

oo 

— 

100 

10 

4.18 

4.10 

4.47 

4.98 

0.04 

0 

100 

2 

4.13 

4.05 

5.27 

4.99 

oo 

— 

2000 

1 

3.98 

4.17 

4.65 

13.10 

oo 

— 

2000 

10 

4.35 

7.94 

4.60 

6.12 
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Table 6 


Spectral indices of F TOT for the flare models of Table 5. Am = m 7 — m x . 


* 2 o 


B.(G) 

&tr / 

m x 

m m 

m 7 

Am 

0.40(b) 

0 

100 

i 

3.21 

3.91 

4.03 

0.83 

0.40(p) 

tt/2 

100 

i 

3.22 

3.91 

4.04 

0.82 

0.40(b) 

0 

2000 

l 

3.21 

3.91 

4.07 

0.86 

0.40(p) 

tt/2 

2000 

l 

3.23 

3.93 

4.32 

1.09 

0.40(b) 

0 

10 4 

l 

3.19 

3.96 

4.65 

1.46 

0.40(p) 

7t/2 

10 4 

i 

3.28 

4.27 

6.00 

2.72 

oo 

— 

100 

l 

3.21 

3.91 

4.03 

0.83 

oo 

— 

100 

2 

3.22 

3.89 

4.03 

0.81 

oo 

— 

100 

5 

3.38 

3.91 

4.06 

0.68 

oo 

— 

100 

10 

3.52 

3.95 

4.11 

0.59 

0.04(b) 

0 

100 

2 

3.19 

3.19 

4.03 

0.84 

oo 

— 

2000 

1 

3.21 

3.91 

4.16 

0.95 

oo 

— 

2000 

10 

3.67 

4.44 

5.01 

1.33 


39 


Figures 


Figure 1.- Electron flux (arbitrary units) vs. pitch angle (degrees) for the collision-only 
case with E = 10.6 MeV, 8 = 5, aj = 0, a\ = 0.40, n e = const. = 10 10 cm -3 , 
and constant B. The analytic results, from equation (15), are given by o’s and the 
numerical results by solid lines. Each curve is labeled by the value of r c jr\. 

Figure 2.- Electron Flux vs. pitch angle at r c = 77, for 16 keV, 1 MeV and 10.6 MeV 

electrons for the collision dominated (R 0 = 3 x 10 -9 ) model with 8 = 5, Q X = 0 and 
= 0.40. Each curve is labeled by the value of the electron energy in MeV. The 
solid lines are the numerical results, and the o’s are the analytic results. The flux 
has been normalized to 1.0 at a = 0. 

Figure 3.- Same as Figure 2, for the pancake injection model (<*1 = 7t/2) with 8 = 5, and 
al = 0.40. There are no analytic results for this case. 

Figure 4.- The Flux integrated over depth F r versus pitch angle for the beam (qi = 0, 
solid lines) and pancake (ai = 7t/2, dashed lines) models with R a = 3 x 10 -9 , 8 = 5 
and ol\ = 0.40, for 300 keV and 10.6 MeV. 

Figure 5.- The flux integrated over pitch angle versus normalized column depth r c /rj 
for the model with a x = 0, a] = 0.40 and R a = 3 x 10 -9 . The lines reperesent 
numerical results and the o’s the analytic results from equation (17). These curves, 
labeled by the value of the electron energy in MeV, are nearly independent of the 
injection parameters such as <5, a*, etc . . . 

Figure 6.- Electron flux versus pitch angle for the synchrotron dominated case with 
R 0 = 3 x 10 3 , at t $ = 1 for an = 0, a* = 0.40, with each curve labeled by 
the value of the electron energy in MeV. Solid lines represent numerical results, 
and the o’s the analytic results from equation (A. 24). The dashed line depicts the 
injected flux. All the fluxes axe normalized to the flux at r, = 0, a = 0. 
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Figure 7.- Flux at zero pitch angle versus energy, for the model with R 0 = 3 x 10 3 . The 
dashed lines are the analytic results of equation (A. 24), and the solid lines are the 
numerical results; all are labeled with the appropriate value of r s . Here the flux is 
multiplied by E 6 , so that the curve for zero depth is flat. 

Figure 8 - The Flux integrated over depth F r versus pitch angle for the beam (solid 
lines) and pancake (dashed lines) models with R a = 1.3, 8 = 5 and a 2 = 0.40, 
and 300 keV and 10.6 MeV electrons. 

Figure 9.- Flux integrated over pitch angle multiplied by E 6 versus electron energy 
for the models with <5 = 5, a 2 0 = 0.40, and R = 1.3. The curves are labeled with the 
appropriate values of t s . The solid lines are for ot \ = 0 (beam) and the dashed lines 
are for O! = 7 t/ 2 (pancake). 

Figure 10.- Electron flux versus pitch angle at the top of the loop ( s = 0) for the ‘short’ 
converging field cases with 8 = 5, a 2 = c© and b = 10. The solid lines axe for the 
case with R 0 ( top) = 3 x 10 -3 at the top, the dashed lines are for R 0 ( top) = 1.3 
and each curve is labeled by the value of the electron energy in MeV. The 16 keV 
results have been shifted by 10 -10 , and the 1 MeV results have been shifted by 
lO" 3 . 

Figure 11- Normailized electron flux versus pitch angle for 10.6 MeV at r c = rj, 8 = 5 
and a 2 = oo. Curves labeled "U" are for the uniform field cases and curves labeled 
"C" axe for the ‘long’ convergence cases with St = 3.7 x 10“ 9 cm -1 . The solid lines 
axe for the models with R 0 = 3 x 10 -9 at the top, and dashed lines are for cases 
with R 0 = 1-3 at the top. 

Figure 12.- Normailzed electron flux versus pitch single for the flare model with 8 = 3.75, 
a 2 = oo, and b = 10. The solid lines are for r c = 0, and the dashed lines are for 
t c = r). Each curve is labeled by the value of the electron energy in MeV. 
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Figure 13.- Electron flux versus pitch angle for 1 MeV electrons for a model with a 

density of Fe ions of 10 27 cm~ 3 and a magnetic field of 10 12 G (or R 0 = 1.2). The 
curves are labeled by the value of r c « t 3 . 
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